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Conventional superconductors are characterized by a single energy scale, the 
superconducting gap, which is proportional to the critical temperature (Tc) (ij 

In hole-doped high-T^ copper oxide superconductors, previous experiments [2 
0] have established the existence of two distinct energy scales for doping levels 
below the optimal one. The origin and significance of these two scales are largely 
unexplained, although they have often been viewed as evidence for two gaps, 
possibly of distinct physical origins {s-Sj. By measuring the temperature de- 
pendence of the electronic Raman response of Bi2Sr2CaCu208+5 (Bi-2212) and 
HgBa2Cu04+5 (Hg-1201) crystals with diff'erent doping levels, we establish that 
these two scales are associated with coherent excitations of the superconducting 
state which disappears at T^. Using a simple model, we show that these two 
scales do not require the existence of two gaps. Rather, a single d-wave super- 
conducting gap with a loss of Bogoliubov quasiparticle spectral weight in the 
antinodal region is shown to reconcile spectroscopic 9-17| and transport 1 18-[21| 
measurements. 

In panels a-b of Fig. 1, we display the Raman spectra x"(a;,T) of Bi-2212 samples with 
different doping levels in Big (antinodal) and i?2g (nodal) geometries (see Methods), for 
several temperatures ranging from well below Tc to 10 K above Tc. In both geometries, these 
spectra show the gradual emergence of a peak as the sample is cooled below T^. An important 
observation is that the Raman responses at Tc and just above Tc are essentially identical. 
This suggests that the emergence of the Big and B2g peaks below Tc is associated with 
coherent excitations of the superconducting state. The Big peak (measured at T = 10 K) 
is seen to decrease in intensity as the doping level is reduced, and disappears altogether at 
the lowest doping p ~ 0.1. In contrast, the i?2gpeak intensity remains sizeable even at low 
doping levels. 

In order to clearly reveal the temperature-dependence of the Big and B2g peaks, we have 
plotted the normalized areas of the -Bi^and -82^ peaks in Fig.le as a function of T/Tc. This 
plot demonstrates that the peak intensities vanish at Tc for all doping levels, providing quan- 
titative support to our interpretation as coherence peaks of the superconducting state. The 
subtracted spectra displayed in Fig. Ic-d (see caption) allow for a quantitative determination 
of the characteristic energy scales associated with the superconducting coherence peaks (see 
plot in Fig. If). It is seen that the Big and B2g energy scales track each other at high doping 
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levels {p ^ 0.19), but depart from each other as doping is reduced. The -Big energy scale 
increases monotonically as doping is reduced, while the energy scale follows a dome-like 
shape approximately similar to that of the critical temperature Tg. 

In order to shed light on the origin of these two energy scales, we consider a very simple 
phenomenological model of a superconductor with a gap function A(0). The angle is 
defined by cos(20) = cos — cos ky and the gap function vanishes at the nodal point A(0 = 
7r/4) = while it is maximal at the antinodes A(</) = 0) = Amax ■ Within a Fermi liquid 
description, ^he quasiparticle contribution to the Raman response in the superconducting 



state reads [3|, : 

^; = ^ (0)z2(0)A2(0)^^M^=\ (1) 

In this expression, ((■ ■ ■))fs denotes an angular average over the Fermi surface, 7^13,^29 
the Raman vertices which read 7^, (0) = 7^ cos 2(f) and 7b2„ (0) = 7b, sin 20, respectively. 



, and A(0)Vv/f^2 _ 4^(0)2 

is a BCS coherence factor. The function Z{(j)) is the spectral 
weight of the Bogoliubov quasiparticles, while A(0) is a Fermi liquid parameter associated 
with the coupling of these quasiparticles to the electromagnetic field. 

In the following, we will show that the angular dependence of the quasiparticle renor- 
malization, Z.A(0), plays a key role in accounting for the experimental observations. In the 
geometry, the Raman vertex 7_Big (0) is peaked at the antinode = 0, resulting in a 
pair-breaking coherence peak at hQsig = 2 Amax due to the singularity of the BCS coher- 
ence factor. The weight of this peak is directly proportional to the antinodal quasiparticle 
renormalization (ZAan)"^ = (ZA)^(0 = 0). Hence, the fact that the coherence peak 
loses intensity at low doping strongly suggests that ZAan decreases rapidly as doping is 



reduced, in qualitative agreement with tunneling |l2l.ll3| and ARPES measurements (221.123 1. 

In the geometry, the situation is more subtle because the Raman vertex is largest at 
the nodes, where the BCS coherence factor vanishes. As a result, the energy of the coherence 
peak depends sensitively on the angular dependence of the quasiparticle renormalization 
ZA(0). If the latter is approximately constant along the Fermi surface, then the energy of 
the -829 peak is determined solely by the angular extension of the Raman vertex 7^3^ (0). 
In contrast, let us consider a ^A(0) which varies significantly, from a larger value ZANsi 
the node to a small value ZA^Arat the antinode, with a characteristic angular extension 07v 
around the node, smaller than the intrinsic width of the Raman vertex 7^29 (0)- Then, it is 
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(pN itself which controls the position of the B2gpeak: h^B2g = 2A(0Ar). As shown below, 
this key feature explains the origin of the differentiation between the two energy scales in 
underdoped cuprates. 

To proceed further in the simplest possible way, we consider a simple crenel-like shape 
for ZA(0), varying rapidly from ZA^r for < < 7r/4 to ZKan < ZK^ for < < 0Ar 



(Fig. 2.AII-CII). Furthermore, we adopt the often-used [2J] parametrization of the gap 
function A(0) = Amaz [i? cos 20+ {1 — B) cos 60], consistent with d-wave symmetry where 
the nodal slope of the gap fA = (9A/90|0=7r/4 = 2(4_B — 3)Am,ax does not necessarily track 

^max ■ 

We thus have 5 parameters: A^ax , the nodal slope of the gap, fAor B (see Methods), 
ZAan , ZAjv and the angular extension 07v • These parameters are determined by attempting 
a semi-quantitative fit to our spectra, obeying the following constraints: 

• The maximum gap A^ax is determined from the measured energy of the Big peak 
according to 2Amax = ^^Si^ 

• The antinodal quasiparticle renormalization ZAan is determined such as to reproduce 
the intensity of the Big coherence peak. 

• The angular extension 0Ar is determined from the energy of the nodal coherence peak. 
Throughout the underdoped regime, this amounts to 2A(0Ar) = hflB2g discussed 
above. 

• The nodal renormalization ZAjyis constrained to insure that the ratio {ZAj^fY/vA 
does not change as a function of doping level, at least in the range 0.1 <]?< 0.16 (as 
shown in Refs. [sl, [lo|). This ratio controls the low-frequency slope of the i?2g Raman 
response. We assume here that the density of states Np (associated with the Fermi 
velocity perpendicular to the Fermi surface) does not depend sensitively on doping 
level in this range. 

These 4 constraints leave one parameter undetermined, which can be taken as the de- 
viation of the gap function from a pure cos — cos ky form, as measured by the ratio 

/ (2 Amax ) = 45 — 3 of the nodal velocity to the gap maximum. We will thus consider 
three possible scenarios: 



4 



• (A) Pure cosA;^^^ — cos fc^ gap: va = ^Amax {B = 1). This corresponds to a super- 
conducting gap involving a single characteristic energy, which increases as the doping 
level is reduced. 

• (B) va tracks the critical temperature Tc. In this case, the gap function is truly 
characterized by two scales varying in opposite manner as the doping level is reduced. 

• (C) f A remains constant as a function of doping. This is also a two-scale supercon- 
ducting gap scenario, although with a milder variation of f a • 

In Fig. 2, we display the Big and i?2g Raman spectra calculated in the framework of this 
simple theoretical analysis, following each of the three scenarios (A-C) above. We observe 
that the main aspects of the experimental spectra, and most importantly the existence of 
two energy scales ^Big , ^B2g varying in opposite manners as a function of doping, can be 
reproduced within any of the three scenarios. 

A common feature between all three scenarios is that the quasiparticle renormalization 
function ZA{(f)) varies significantly along the Fermi surface. Quasiparticles have a large spec- 
tral weight ZA/vronly on a restricted region around the nodes, defined by (p^, corresponding 
to a fraction fc = (vr/4 — 0Ar)/(-7r/4) of the Fermi surface. While A^ax increases with falling 
doping, A((/)7v) decreases because of the rapid contraction of the 'coherent fraction' fc , lead- 
ing to the opposite doping dependence of the two scales, as illustrated on Fig. 3a. We note 
that linearizing the gap function in the coherent region is a reasonable approximation for 
the A and C scenarios, leading to the relation hQB2g = f /c'^a ksTc which links the nodal 
{B2g) energy scale (proportional to Tc), the nodal velocity and the coherent fraction. This 
approximation is not valid for scenario B, because linearization is not accurate there. 

It is clear that having uniformly coherent Bogoliubov quasiparticles along the Fermi 
surface is inconsistent with our data, especially in view of the rapid suppression of the 
Big coherence peak and the corresponding decrease of ZKan ■ This is one of the key conclu- 
sions of the present work. 

In panels (A.V-CV) of Fig. 2, we display a calculation of the corresponding tunneling 
conductance, (see Methods). One sees that the two energy scales have clear signatures in 



the tunneling spectra, the nodal (i?2g ) one correspondin 
a feature has indeed been observed in tunneling spectra 



to a 'kink'-like feature. Such 
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17|, although it has been 



usually interpreted as evidence that and A^ax are two distinct scales (scenarios B, or 



C), without considering the effect of Z{(j)). Our work demonstrates that the STM 'kink' 
results quite generally from the rapid decrease of the quasiparticle spectral weight as one 
moves from nodes to antinodes. 

Although the above features are common to all three scenarios, there are two key differ- 
ences between them. The first one is qualitative: in scenario A, ZAj^j increases as doping 
level is reduced, while it decreases for scenario B and stays constant for scenario C. The sec- 
ond, quantitative, difference is the rate at which the coherent fraction of the Fermi surface 
/c decreases with underdoping, being largest for scenario A and smallest for B (Fig. 3b). 

A determination of the coherent fraction /c^*" has been reported from heat-capacity 
HC) measurements fisl . 1^, as reproduced on Fig. 3b. It was also reported from ARPES 
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25l | in the normal state, that the Fermi arcs shrink upon cooling as ~ T/T*. The doping 
evolution of the coherent fraction fc ^^^^^ oc T^/T* at Tc is displayed in Fig. 3b. Remarkably, 
we find that there is a good quantitative agreement between the doping dependence of 
/c reported from HC and ARPES and our determination from Raman within scenario A 
(a single gap scale v/^ oc A^ax), which thus appears to be favored by this comparison. 
Although this quantitative agreement should perhaps not be overemphasized in view of 
the uncertainties associated with each of the experimental probes, we conclude that this 
single-gap scenario (A) stands out as the most likely possibility. Our interpretation also 
reconciles the distinct doping dependence of the two energy scales with thermal conductivity 



measurements in underdoped samp 
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21| (see, however 



es which, interpreted within the clean limit, suggest that 



26|). 



As noted before, with a single superconducting gap, the relation between the critical 
temperature (or ^B2g ) the coherent fraction reads: ksT^ oc fc^max ■ We are then led 
to conclude that it is the suppressed coherence of the quasiparticles that sets the value of 
Tc, while Amax increases with underdoping. From our findings, we anticipate that the loss 
of coherent of Bogoliubov quasiparticles on a restricted part of the Fermi surface only is a 
general feature of how superconductivity emerges as holes are doped into a Mott insulating 
state. 
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METHODS 



Details of the experimental procedure 

The Bi-2212 and Hg-1201 single crystals have been grown by using a floating zone method 
and flux technique respectively. The detailed procedures of the crystal growth are described 



elsewhere 
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al from Ref. 



|. The doping value p is inferred from using equation of Presland et 
29| : 1 — Tc/TJ""^ = 82.6 (p — 0.16)^. Tc has been determined from magnetic 
susceptibility measurements for each doping level. Raman experiments have been carried 
out using a triple grating spectrometer (JY-T64000) equipped with a nitrogen cooled CCD 
detector. All the measurements have been corrected for the Bose factor and the instrumental 
spectral response. The Big and geometries have been obtained using crossed light 
polarizations at 45° from the Cu — O bond directions and along them respectively. In these 
geometries we probe respectively, the antinodal (AN) and nodal (N) regions correponding 
to the principal axes and the diagonal of the Brillouin zone. 

Calculation of the tunneling conductance 

^ oc Np Z{(j)) (^^Y-''irY-A{(t,)'^ ) where F = ^max /30 is a small damping factor, and 
we have assumed for simplicity that Z{(j)) and ZK{(f)) can be identified. 

ACKNOWLEDGEMENTS 

We are grateful to J. Mesot, A. Yazdani, A. J. Millis, P. Monod and C. Ciuti for very 
helpful discussions. Correspondance and requests for materials should be addressed to A.S. 
(alain.sacuto@univ-paris-diderot.fr). 



[1] Bardeen, J., Copper, L.N, & Schrieffer, J.R. Theory of superconductivity. Phys.Rev. 106, 162 
(1957). 

[2] Deutsche!, G. et al. Coherence and single-particle excitations in the high-temperature super- 
conductors. Nature 397,410 (1999). 



7 



[3] Le Tacon, M. et al. Two energy scales and two distinct quasiparticle dynamics in the super- 
conducting state of underdoped cuprate. Nature Phys. 2, 537 (2006). 

[4] Tanaka, K. et al. Distinct Fermi-momentum-dependent energy gaps in deeply underdoped 
Bi2212. Science 314, 1910 (2006). 

[5] Lee, W. S. et al. Abrupt onset of a second energy gap at the superconducting transition of 
underdoped Bi-2212. Nature 450, 81 (2007). 

[6] Boyer, M. C. et al. Imaging the two gaps of the high-temperature superconductor 
Bi2Sr2CuOfi+^. Nature Phys. 3, 802 (2007). 

[7] Millis, A. Gaps and Our Understanding. Science314, 1888 (2006). 

[8] Hufner, S., Hossain, M. A., Damascelli, A. & Sawatzky, G. A. Two gaps make a high- 
temperature superconductor. Rep. Prog. Phys. 71, 062501 (2008). 
[9] Guyard, W.et al. Breakpoint in the evolution of the gap through the cuprate phase diagram. 
Phys. Rev. B 77, 024524 (2008). 
[10] Blanc, S. et al. Quantitative Raman measurement of the evolution of the Cooper-pair density 

with doping in Bi2Sr2CaCu20s+5 superconductors. Phys. Rev. B 80, 140502R (2009). 
[11] Devereaux, T. P. & Hackl, R. Inelastic light scattering from correlated electrons Rev. Mod. 
Phys. 79, 175 (2007). 

[12] McElroy, K.et al. Coincidence of Checkerboard Charge Order and Antinodal State Deco- 
herence in Strongly Underdoped Superconducting Bi2Sr2CaCu20^j^s. Phys. Rev. Lett. 94, 
197005 (2005). 

[13] Kohsaka, Y. et al. How Cooper pairs vanish approaching the Mott insulator in 
Bi2Sr2CaCu208+5- Nature 454, 28 (2008). 

[14] Kanigel, A. et al. Protected Nodes and the Collapse of Fermi Arcs in High-Tc Cuprate Su- 
perconductors. Phys. Rev. Lett. 99, 157001 (2007). 

[15] Kondo, K., Khasanov, R., Takeuchi, T., Schmalian, J. &: Kaminski, A. Competition between 
the pseudogap and superconductivity in the high-Tc copper oxides. Nature 457, 296 (2009). 

[16] Chatterjee, U. et al. Observation of a d-wave nodal liquid in highly underdoped 
Bi2Sr2CaCu20s+5. Nature Phys. 5, 1456 (2009). 

[17] Pushp, A. et al. Extending Universal Nodal Excitations Optimizes Superconductivity in 
Bi2Sr2CaCu208+5. Science 26, 1689 (2009). 

[18] Wen, H.H. et al. Pseudogap, superconducting energy scale, and Fermi arcs of underdoped 

8 



cuprate superconductors. Phys. Rev.B 72, 134507 (2005). 
[19] Matsuzaki, T., Momono, N., Oda, M. & Ido,M. Electronic Specific Heat of La2-xSrxCu04: 

Pseudogap Formation and Reduction of the Superconducting Condensation Energy J. Phys. 

Soc. Jpn. 73, 2232 (2004). 
[20] Sutherland, M.et al. Thermal conductivity across the phase diagram of cuprates: Low-energy 

quasiparticles and doping dependence of the superconducting gap. Phys. Rev.B 67, 174520, 

(2003). 

[21] Hawthorn, D.G et al. Doping dependence of the superconducting gap in Tl2Ba2CuOQ-^.s from 
heat transport. P/iys. Rev. B 75, 104518 (2007). 

[22] Ding, H. et al. Coherent Quasiparticle Weight and Its Connection to High- Tc Superconduc- 
tivity from Angle- Resolved Photoemission, Phys. Rev. Lett. 87, 227001 (2001). 

[23] Vishik et al. A momentum-dependent perspective on quasiparticle interference in 
Bi2Sr2CaCu20s+5 Nature phys. 5, 718, (2009). 

[24] Mesot, J. et al. Superconducting Gap Anisotropy and Quasiparticle Interactions: A Doping 
Dependent Photoemission Study. Phys. Rev. Lett. 83, 840, (1999). 

[25] Kanigel, A. et al. Evolution of the pseudogap from Fermi arcs to the nodal liquid. Nature 
Phys. 2, 447 (2006). 

[26] Sun, X. F. et al. Electronic Inhomogeneity and Breakdown of the Universal Thermal Conduc- 
tivity of Cuprate Superconductors, Phys. Rev. Lett. 96, 017008, (2006). 

[27] Wen, J. S. et al. Large Bi-2212 single crystal growth by the floating-zone technique. J. of 
Crystal Growth. 310, 1401 (2008). 

[28] Bertinotti A. et al. In studies of High Temperature Superconductors Vol.23 (ed. A. Narlikar), 
27-85 (Nova Science Pubhshers New York, (1997). 

[29] Presland M. R. et al. General trends in oxygen stoichiometry effects on Tc in Bi and Tl 
superconductors. Physica C 176, 95 (1991). 



9 



3^ 



- njf 10K 75K 

jOfy 40K 85K 

60K 


b 

p 


=0.21 Tj,= 75K 


10K ^ 

40K f\ 

70K 

84K r 


// 
^ P 


=0.19 T=84K 


. 10K 

40K A. 

70K 

91K f 

100K A^~"'""v^'!W 




1 

=0.16 T^=91K 


. 10K 

4QK 

6QK 

87K / \ 

95K /''"NSw 

^^^^"^ 


K 


=0.14 Tj,= 87K 


. 10K 

40K 

60K 

75K 

85K /\ 




=0.12 Tj.= 75K 


10K 

40K 

55K 

68K 

78K 


/ p 


=0.1 Tj,= 68K 



g 



3 



200 400 600 800 

CO (cm') 



200 400 600 800 

CO (cm') 





d A B^g 




" : p=o.2i 




A 




^ : p=o.i9 






V : p=o.i6 
















^^^^^^^ 
: p=o.io 




) 200400 600 800 


200 400 600 800 


CO (cm') 


CO (cm') 



Si.o 



CD 
D. 



T3 

CD 
N 



0.6 - 



CO 0.4 

E 
o 

2 0.2 - 



0.0 



Bi-2212 B^g 

■ 0.21 A 0.21 

■ 0.19 A 0.19 

■ 0.16 A 0.16 
0.14 0.14 

■^fl a" ^ ■ ° 

A ^ ^ A 0.10 



1^ A 

? a' 



1.0 \ fl A 

Hg-1201 ^ 



0.5 



■22 A 



1.0 



■ 4 



0.0 0.2 0.4 O.f 

T/T 



0.8 1.0 




0.08 0.12 0.16 0.20 0.24 

Doping p 



FIG. 1: (color online) Temperature dependence of the Raman spectra. a,b: Raman spectra, 
x"{to,T) of Bi-2212 single crystals for several doping levels in Big (Antinodal) and B2g (Nodal) 
geometries, c, d: Raman spectra subtracted from the one measured at 10 K above Tc for each 
sample in each geometry. A direct visual comparison of the subtracted spectra to a reference energy 
(chosen as the peak position for the most overdoped sample, drawn as a guide to the eyes) clearly 
reveals the distinct doping-dependence of these two energy scales, e: Temperature dependence of 
the normalized areas of the Big and B2g peaks with respect to the area measured at T = 10 K. 
The inset displays these areas for Hg-1201 crystals, f: Energy scales associated to the Big and i?2g 
peaks for both Bi-2212 and Hg-1201 crystals. The dashed lines are guides to the eyes and track 
the locations of the superconducting peak maxima. 
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FIG. 2: (color online) Three scenarios for underdoped cuprates. A.I-C.I: Doping evolution 
of the superconducting gap in the three scenarios (A-C). A.II-C.II: Angular dependences of the 
quasiparticle spectral weights ZA(0) as a function of doping level for each scenario. The angular 
dependence of the B2g Raman vertex is shown in dotted line (see All). ATII-C.III, ATV-C.IV: 

Calculated Raman spectra for each scenario in Big and B2g geometries. The locations of the Big 
and B2g peaks are respectively controlled by the gap energy at the antinodes <I> = and at the 
angle (^n- Note that the low energy slope of the i?2g Raman response is constant with doping as 



observed experimentally 



10l |. A.V-C.V: Tunneling spectra for the three scenarios. 
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FIG. 3: (color online) Fermi surface coherent fraction and a single d-wave superconduct- 
ing gap. a: Scenario A where coherent Bogoliubov quasiparticles are partially suppressed on 
restricted parts of the Fermi surface and the superconducting gap has a single d-wave shape, b: 
Doping evolution of the coherent fraction of the Fermi surface (/c) in the three scenarios (A-C). 



The fc^^^^^ curve is deduced from Larc/Lfuii{Tc) = 1 — 0.70^ (see Ref |25l|). expressed as a 
function of the doping level, fc is extracted from the zero temperature specific heat coefficient 



7(0)n: 
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